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ABSTRACT 

The millimeter wavelength emission from GRB 991208 is the second brightest 
ever detected, yielding a unique data set. We present here well-sampled spectra 
and light curves over more than two decades in frequency for a two-week period. 
This data set has allowed us for the first time to trace the evolution of the 
characteristic synchrotron self-absorption frequency u a and peak frequency u m , 
and the peak flux density F m : we obtain u a oc £-° 15 ± 012 ; Vm ^ t -L7±0 - 4 , and 
F m oc t _0 - 47±011 . From the radio data we find that models of homogeneous or 
wind-generated ambient media with a spherically symmetric outflow can be 
ruled out. A model in which the relativistic outflow is collimated (a jet) can 
account for the observed evolution of the synchrotron parameters, the rapid 
decay at optical wavelengths, and the observed radio to optical spectral flux 
distributions that we present here, provided that the jet transition has not 
been fully completed in the first two weeks after the event. These observations 
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provide additional evidence that rapidly decaying optical/X-ray afterglows are 
due to jets and that such transitions either develop very slowly or perhaps never 
reach the predicted asymptotic decay F(t) oct~ p . 



Subject headings: gamma rays: bursts - radio continuum: general - cosmology: 
observations 



1. Introduction 

GRB 991208 was detected with the Interplanetary Network (IPN) by the Ulysses, 
WIND and NEAR spacecraft on December 8, 1999, at 04:36:52 UT. The gamma-ray 
properties, the localization and the subsequent discovery of the radio counterpart to 
GRB 991208 are presented by Hurley et al. (2000). 

The IPN can typically localize events within one or two days. Therefore, important 
spectral transitions, such as the passage of the synchrotron peak and cooling frequencies 
( [Ramaprakash et at. 199~8| , |Galama et at. 1998| ), that occur at optical and infrared 



wavelengths a few hours after the event are missed by optical and infrared follow-up 
observations on IPN bursts. Fortunately, the radio range allows similar studies on time 
scales from days to weeks. 

Here we present an extensive set of radio observations monitoring GRB 991208 between 
1.4 and 250 GHz during the first two weeks after the event. It is the second brightest 
afterglow detected at millimeter wavelengths to date ( [Shepherd et al. 1999 ). The brightness 



of the afterglow allows us to trace the evolution of the characteristic synchrotron peak 
frequency, the self-absorption frequency, and the peak flux density. The observations and 
results are presented in §|2], and compared with appropriate models in §|3|. 



2. Observations and results 

Observations were made from 1.43 to 250 GHz, at a number of facilities: the Max- 
Planck Millimeter Bolometer (MAMBO; [Kreysa et al. 1998| ) at the IRAM 30-m telescope 



on Pico Veleta (Spain; [Baars et al. 1987| ), the IRAM Plateau de Bure Interferometer (PdBI) 
in the French Alps ( |Guilloteau et al. 1992| ), the Owens Valley Radio Observatory (OVRO) 
millimeter-wave array, the OVRO 40-meter telescope, the Ryle Telescope at Cambridge 
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(UK), and the NRAO Very Large Array (VLA)P]. We have detailed our observing and 
calibration strategy in Kulkarni et al. (1999), Galama et al. (1999), Frail et al. (2000), and 
Berger et al. (2000). A log of the observations is provided in Table [I]. 



The millimeter emission from GRB 991208 is one of the brightest ever detected (|Bremer 



et al. 1998 , |Smith et al. 19991 , [Shepherd et al. 1998| ) enabling us to obtain a unique data 



set. Rather than single-epoch snapshot spectra (|Galama et al. 1998|) , we have well-sampled 
spectra and light curves over more than two decades in frequency for a two-week period. 
GRB 991208 is therefore well-suited to a study of the broad-band evolution of the radio 
afterglow in the first two weeks after the burst. We display light curves for a subset of these 
data in Fig. [TJ. The data from Table [TJ at 22 GHz, 30 GHz and 86 GHz are too sparse to 
plot in this form. 

Several trends are apparent. At the highest frequencies (250 GHz, 100 GHz and 15 
GHz) the flux densities are either declining or constant during the first week, but thereafter 
fade below detectability. At frequencies below 10 GHz we see erratic flux density variations 
due to interstellar scintillation (ISS; see e.g. |Goodman 1997| ; details on the ISS properties 
of GRB 991208 can be found in Galama et al. 2000). Therefore, below 10 GHz it is harder 
to discern a pattern. The overall trend in Fig. [I] is for the peak flux density to decline with 
decreasing frequency, while the time-to-maximum increases. 

An alternate way to view the data in Table [I] is to construct instantaneous spectral 
flux distributions at several epochs, which were chosen to obtain a maximum number 
of observations sufficiently close in time. All values were brought to the same epoch by 
applying a correction that was determined by fitting cubic splines to each of the light curves 
in Fig. [TJ. We have added in quadrature an additional error of 25 % of the synchrotron 
spectral model flux (see §2.1] ) to the fluxes below 10 GHz to reflect the uncertainty due to 
ISS (see |Galama et al. 2000] ) . The resulting spectral flux distributions are presented in Fig. 



[| The spectra for all four epochs show the same overall morphology - flat from 10 GHz to 
250 GHz, and dropping precipitously below 10 GHz. 

In addition, we reconstruct the radio to optical spectral flux distribution at the epoch 
of December 15.5 UT by including K, R and I-band detections ( [Bloom et al. 1999] , [Masetti| 
et al. 199"9"|, |Sagar et al. 2000]) . All values were brought to the same epoch by applying a 



correction using the slope of the fitted optical light curve ( [Sagar et al. 2000] ). We corrected 



for Galactic foreground reddening (A_r = 0.043, Aj = 0.031, and A^ = 0.006; as inferred 



11 The NRAO is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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from the dust maps of Schlegel, Finkbeiner and Davis 19980). The result is presented in 
Fig. 0. 



2.1. Spectral fits 

The cm/mm/optical to X-ray afterglow emission is believed to arise from the forward 
shock of a relativistic blast wave that propagates into the circumburst medium (see [Piran 



1999, van Paradijs, Kouvcliotou fc Wijcrs 2000 for reviews). The synchrotron afterglow 



spectrum has four distinct regions: F v oc v 2 below the synchrotron self-absorption frequency 
u a (y < z/ a ); F u oc z/ 1//3 up to the synchrotron peak at u m (v a < v < u m ); F v oc i>~(p~ 1 ^/ 2 
above the peak up to the cooling frequency u c (u m < v < u c ); and F v oc v~ p l 2 above cooling 
[y > v c \ see e.g. pari, Piran fc Narayan 1998| ; we have assumed slow cooling, i.e. v m < z/ c ). 
Here, p is the power-law index of the electron energy distribution. 

The spectral flux distribution of Dec 15.5 UT (Fig. |3|) peaks at z/ m ~ 40 GHz. A fit 
to the radio data is provided by a synchrotron spectrum from a relativistic blast wave as 
specified by Granot, Piran and Sari (1999a; see their Fig. 10 for the equipartition field 
model). We scale their dimensionless functional form by peak frequency h> m and peak flux 
density F m to derive a function g(u) with asymptotic behavior of v 1 ^ {y « z/ m ) and 
^-(p- 1 )/ 2 >> j, m )_ \y e account for synchrotron self-absorption at v & by multiplying g{y) 
by F v — [1 — exp(— t)]/t, where r = (^/z/ a ) -5 / 3 QGranot, Piran fc Sari 1999b|) . 



First we fit the radio to optical spectral flux distribution of December 15.5 UT and 
find: p = 2.524 ± 0.016, z/ a = 3.6±J;| GHz, u m = 35 ± 4 GHz, and F m = 2.85 ± 0.15 mJy 
(Xi — 1-6; 6 degrees of freedom; d.o.f.). The spectrum is well described by synchrotron 
emission from a power-law electron distribution with index p = 2.52, self-absorbed at low 
frequencies; the fit is shown in Fig. 0. 

Next, we fix p = 2.52 and fit all the epochs of the radio spectral flux distributions (the 
fits are shown in Fig. 0): the derived values of u a , v m and F m and their uncertainties are 
given in Table |2|. These are plotted in Fig. ^, where a clear temporal evolution of u m and 
F m is apparent. We note that assuming p = 2.2 (see §0) does not change the results of 
the fits (i.e., when no optical data is included the fits are rather insensitive to the assumed 
value of p). To characterize this evolution, we made power-law least squares fits to the data. 
The results are given in Table 0. 



sec http: / /astro. berkeley.edu/davis / dust /index. htm] 



- 5 - 



3. Discussion 

The brightness of the radio counterpart to GRB 991208 has allowed for the first time 
to trace the evolution of the characteristic synchrotron parameters z/ a , v m and F m . We may 
now compare this result directly with expectations from relativistic blast wave models. 

Currently popular scenarios for the origin of GRBs are: a compact object merger (the 
neutron star-neutron star, [Eichler et al. 198~9|, and neutron star-black hole merger models, 



Lattimer fc Schramm 1571 , |JNarayan, Paczyhski fc Fir an 19933 ), and the core collapse of 



a very massive star ('failed' supernova or hypernova |Wooslcy 1993| , Paczyhski 1998 ). In 



principle afterglow observations can distinguish a compact object merger, which is expected 
to occur in a constant density interstellar medium (constant density ISM model; see |Sari,| 
Piran fc Narayan 1998| and [Wijers fc Galama 1999| for details), from a hypernova, where 



the circumburst environment will have been influenced by the strong wind of the massive 
progenitor star (wind model; see |Uhevalier fc Li 1999a| , Chevalier fc Li 1999b| ). Below we 



compare each of the afterglow models separately with the observations. The expected 
scalings (we assume adiabatic expansion of the remnant) for each of the models, including 
the case where the outflow is collimated in a jet, are also provided in Table |3| for comparison 
with the evolution derived from the spectral fits. 

The observed decay of the peak-flux density F m is inconsistent with the constant 
density ISM model (see Table |J), which predicts a constant peak flux. Also, the optical 
afterglow of GRB 991208 (R band) displays a rapid decay: a = —2.2 ± 0.1, F u oc t a flSagar 



et al. 2000] ; see also |Castro-Tirado et al. 2000|) . A constant density ISM model would 



require p > 3.6 in order to account for the rapid decay observed at optical wavelengths. For 
p > 3.6 the synchrotron spectrum above the peak should be much steeper than observed 
(F u oc z/~( p_1 )/ 2 ). We note that optical extinction local to the host would cause the observed 
spectrum to be even shallower. The constant density ISM model can thus be ruled out 
on the basis of the spectral flux distribution on 1999 Dec. 15.5 UT, which shows a flat 
spectrum corresponding to p — 2.52 (see Fig. ^). 

Comparing GRB 991208 with the wind model (see Table |]) we find that the evolution 
of the peak-flux density F m and the peak frequency expected, but we find no 

evidence for a decrease with time of the self-absorption frequency z/ a . And, similar to the 
case of the constant density ISM model, a wind model would require p > 3.3 in order to 
account for the rapid decay observed at optical wavelengths. We therefore also reject the 
wind model. 



Finally, the outflow may be collimated into jets (jet model; see for details [Rhoads 



1999| , [Sari, Piran fc Halpern 1999| , Panaitescu fc Meszaros 1999|) . Once sideways expansion 



- 6- 



of the jet dominates the dynamics of the relativistic-blast wave, a jet plus constant ambient 
density medium (ISM + jet) or a jet plus wind stratified ambient medium (wind + jet) will 
have very similar properties. In this stage, the morphology of the ambient medium will be 
difficult to constrain. 

The observed evolution of u a , and Z^ m £1X6 clS expected in the jet model, but the drop in 
peak-flux density F m is not as rapid as predicted (see Table ||). To obtain the predicted 
decay rate F m oc t^ 1 we have to raise x 2 by 21.0 (corresponding to > 4cr); the observed 
decay is thus significantly different from the predicted decay. And, in the case of a jet, at 
late times, when the evolution is dominated by the spreading of the jet, the model predicts 
a = —p, and the R-band measurement thus corresponds to p = 2.2 ± 0.1. For p = 2.2 the 
synchrotron spectrum above the peak should be shallower than observed. But, we find that 
p = 2.52 gives a satisfactory description of the Dec 15.5 UT data (Fig. |3p. For example 
assuming p = 2.2 provides a worse fit (Xr = 3-4 as opposed to 1.6 for p = 2.52; Dec. 15.5 
UT radio to optical spectral flux distribution; Fig. |3|). 

A p = 2.2 jet model could only be saved if one assumes a cooling break between the 
mm and optical passbands, around v c = 5 x 10 13 Hz (the optical decay would still satisfy 
a = —p, i.e. p = 2.2 ±0.1). However, in that case the optical/infrared synchrotron spectrum 
F^ oc v~ v l 2 = v" 1 - 1 , i.e. significantly steeper than observed {(3 = —0.75 ± 0.03, F u ex u 13 ; 
Sagar et al. 2000]) . We therefore reject the possibility that p = 2.2. 



In the jet model, the shallow decay of the peak-flux density F m and the apparent 
discrepancy in the value of p derived from the spectrum and from the light curve, can 
be understood if during the first two weeks the GRB remnant was still undergoing 
the transition from a quasi-spherical (either constant density ISM or a wind), i.e. 
F(t) oc r( 3 P- 3 )/ 4 ~ r L1 (constant density ISM), or F(t) oc H 3 ^ 1 )/ 4 ~ r 1 - 6 (wind) to a jet 
evolution, F(t) oc t~ p ~ t~ 2 b . As discussed by Kumar and Panaitescu (2000) this transition 
takes place over at least an order of magnitude, but more likely several orders of magnitudes, 
in observer time. Since typically we expect the jet transition to take place around one to 
several days (pari, Piran fc Halpern 1999|) it is natural to expect that the transition was still 



ongoing; the evolution of the synchrotron parameters will then be somewhere between fully 
spherical at early times (either constant density ISM or a wind) and fully dominated by a 
jet later. This is consistent with the data, would explain the shallow decay of the peak flux 
density, and would also explain why the optical decay rate exponent a = —2.2 ( [Sagar etal. 



|2000[ ) is not as steep as for a fully developed jet transition a = — p = —2.52. Determination 
of p from light curves requires observations after the transition has been completed and 
from the optical data we cannot establish whether this is the case. Also the light curve is 
not necessarily expected to approach the asymptotic limit a = —p. The value of p derived 
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from the optical light curve is thus rather uncertain and the value of p = 2.52 that we have 
inferred from the spectral flux distribution is the more robust number, since p so derived is 
independent of the hydrodynamical evolution of the blastwave. 

A similarly slowly developing jet transition was observed for GRB 990123 ( [Kulkarni ej 
al. 1999a] , [Fruchter et al. 1999| , |Castro-Tirado et al. 1999| ), where after the jet break (at 



tj ~ 2 days) the optical decay rate exponent a ~ —1.8; whereas a steep a = —p = — 2.4 
would have been expected for a fully developed jet transition. Also for GRB 991216 a slow 
transition is seen, but due to host contamination, the late-time decay rate exponent can 
only be constrained to lie between —2.1 and —1.5 ( [Halpern et al. 2000| ). On the other hand 
GRB 990510 and GRB000301C appear to be examples of relatively fast and fully developed 
jet transitions ( jHarrison et al. 1999| , |Berger et al. 2000] ; but see [Kumar fc Panaitescu 2 



In conclusion, the jet model can account for the observed evolution of the synchrotron 
parameters of GRB 991208, the rapid decay at optical wavelengths and for the observed 
radio to optical spectral flux distributions, provided that the jet transition has not been 
fully completed in the first two weeks after the event. These observations provide additional 
evidence that rapidly decaying afterglows are due to jets ( [Sari, Piran fc Halpern 1999]) and 



that at least in some afterglows these transitions are not fully developed. 

Around 10 days after the event u m passed z/ a (u m < v < z/ a ; see Fig. f|) at ~ 8 
GHz; from then on z/ a is expected to evolve differently: v a oc t~ 2 ( 1+ p)/( 4 +p) = £- 14 / 13 (f or 
p = 2.5); we have assumed a jet-model evolution. Also after u m has passed 8.46 (~ 12 days) 
and 4.86 GHz (~ 17 days), the flux densities at these frequencies should decay rapidly 

oc t~ 2,2 ~ 2 5 . Continued monitoring with the VLA flGalama et al. 2000| ) will allow tests of 
these expectations. 



Research at the Owens Valley Radio Observatory is supported by the National Science 
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Sargent for generously allocating the time necessary to track the millimeter evolution. We 
thank A. Diercks for useful discussions. 
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Table 1. Radio Observations of GRB 991208 



Date 


Freq. 


Flux 


Remai 


1999 Dec. (UT) 


(GHz) 


(A*Jy) 




11.77 


250 


2600 ± 800 


MAMBO, r atm 


12.38 


250 


2000 ± 500 


MAMBO, r atrn 


14.35 


250 


2000 ± 700 


MAMBO, r atm 


21.38 


250 


100 ± 600 


MAMBO, r atrn 


11.63 


100 


3600 ± 800 


OVRO MA 


13.84 


100 


3500 ± 800 


OVRO MA 


15.72 


100 


1900 ± 800 


OVRO MA 


16.89 


100 


1700 ± 800 


OVRO MA 


15.62 


86.24 


2500 ± 500 


PdBI 


13.63 


30.0 


2700 ± 1300 


OVRO 40-m 


14.62 


30.0 


1700 ± 1000 


OVRO 40-m 


14.84 


30.0 


3600 ± 1200 


OVRO 40-m 


11.51 


15.0 


2100 ± 300 


Ryle 


12.44 


15.0 


2100 ± 400 


Ryle 


13.55 


15.0 


3100 ± 600 


Ryle 


15.61 


15.0 


2300 ± 400 


Ryle 


20.68 


15.0 


300 ± 600 


Ryle 


23.60 


15.0 


500 ± 500 


Ryle 


21.96 


22.49 


1100 ± 300 


VLA 


21.96 


14.97 


1650 ± 250 


VLA 


10.92 


8.46 


707 ± 39 


VLA 


12.04 


8.46 


1630 ± 51 


VLA 


13.81 


8.46 


1900 ± 66 


VLA 


15.91 


8.46 


1990 ± 33 


VLA 


17.81 


8.46 


999 ± 50 


VLA 


18.81 


8.46 


1182 ± 52 


VLA 


19.92 


8.46 


1930 ± 53 


VLA 


20.92 


8.46 


1173 ± 49 


VLA 


21.96 


8.46 


1146 ± 66 


VLA 


10.92 


4.86 


327 ± 45 


VLA 


12.04 


4.86 


622 ± 66 


VLA 


21.96 


4.86 


820 ± 52 


VLA 


15.90 


1.43 


250 ± 60 


VLA 


16.80 


1.43 


250 ± 70 


VLA 


21.96 


1.43 


18 ± 56 


VLA 



0.16 
0.22 

0.54-0.77 
0.16 



- 11 - 



Table 2. Results of Spectral Fitting for GRB 991208 : the epoch of the spectral flux 
distribution, the synchrotron self-absorption frequency u a , the peak frequency u m , the peak 
flux density F m , the reduced chisquared \i an d the degrees of freedom (d.o.f.) of the fit. 



Date 






F m 


X? 


d.o.f 


(UT) 


(GHz) 


(GHz) 


(mJy) 






1999 Dec 11.50 


n.o3:§ 


55tH 


3.9 ±0.4 


0.0 


2 


1999 Dec 13.50 


7±t 


35i| 


3.9 ±0.5 


0.8 


1 


1999 Dec 15.50 


o 0+1.2 


30t* 2 


2.80 ±0.30 


1.0 


3 


1999 Dec 21.50 


8.8i};° 


4.8±?;1 


2.10 ±0.25 


1.2 


3 



Note. - (1) In the spectral fitting we have included the non detections at 1.4 and 250 
GHz on 1999 Dec 21.5 UT as the peak flux density at the location of the afterglow. (2) For 
Dec 21.5 UT we have also fitted using r = 13/6 (see § |2.1| ) to account for the fact that the 
peak frequency v m is below the self- absorption frequency v & at this epoch and so F^ cx z/ 5//2 . 
The result of the fit was the same within the error bars. 
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x 
2 



o.i 



2 



o.i 



x 
2 



o.i 







r^^^^^ 240 GHz : 

\ 1 1 1 1 1 1 1 : 


: — — — j^^^^^^^^^ 

t 1 1 — i — i — i — i — i — ; 


: 15 GHz : 

= : 


: 8.5 GHz : 

• • • • 




\ 1 1 1 1 1 1 1 : 


= — — •* 

• 

b 1 1 1 1 — 1 — 1 — 1 — ■ 


: 4.86 GHz : 

r • 

- ■ : 


: 1.43 GHz : 

r T i 


: 1 



5 10 2 

time (clays) 



5 10 
time (days) 



Fig. 1. — Radio, mm and sub-mm light curves of GRB 981208 from 1.4 to 250 GHz. These 
data are compiled in Table [I] and are taken at various telescopes which include the IRAM 
30-m (250 GHz), the OVRO Millimeter Array (100 GHz), the Ryle Telescope (15 GHz), and 
the VLA (8.5, 4.9 and 1.4 GHz). At 15 GHz there is one additional point on the plot (open 
circle) that was taken with the VLA. All upper limits are plotted as the peak flux density 
at the location of the afterglow plus two times the rms noise in the image. Frequencies 
below 10 GHz show erratic flux density variations due to ISS. Shown are the lightcurves that 
result from assuming an evolution of z/ a , u m , and F m as derived in §0, i.e. u a oc £- - 15±0 - 12 ; 
v m oc r 1,7±0 - 4 , and F m oc t-° A7±0 - u and using the Granot et al. (1999a, 1999b) formulation. 
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Table 3. Inferred evolution of the synchrotron parameters z/ a , u m , and F m for GRB 991208 
and the predicted scalings from the constant density ISM, wind and jet model. 



/A, 



Observed 12.0±l;f xl0^ b s 15±a12 Hz Q.ltH xlO 11 t£ 7±DA Hz 7.2±i:| *-£ 47±CU1 mJy 

( Xr 2 = 5.5; 2 d.o.f.) (x 2 = 1.6; 2 d.o.f.) ( X 2 = 1.2; 2 d.o.f.) 

ISM octl oc^ oc^l 

wind octj^ 5 «t ob 3 / 2 oct^ 1 / 2 

jet oc t obs oc t obs cx t obs 



Note. — (1) t obs is the observer time in days after the event. (2) The scalings for the jet 
model are valid if the evolution is dominated by the sideways expansion of the jet (i.e. at 
late times). 
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Fig. 2. — The radio spectral flux distribution of GRB 991208 from 1.4 to 250 GHz at several 
epochs. Shown are fits of the synchrotron spectra to the radio data from a relativistic blast 
wave as specified by Granot et al. (1999a, 1999b) for p = 2.52. 
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GRB 991208 Dec 15.5 UT 




freq (Hz) 

Fig. 3.— The radio to optical spectral flux distribution of GRB 991208 at Dec 15.5 UT 
(7.3 days after the event). Shown is a fit of the synchrotron spectrum to the radio and 
infrared/optical data from a relativistic blast wave as specified by Granot et al. (1999a, 
1999b) for p = 2.52. 
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Fig. 4. — Evolution of the the synchrotron self-absorption frequency the synchrotron 
peak frequency u m and the peak flux density F m as derived from broad band radio (from 1.4 
to 250 GHz) spectral fits of GRB 981208 at several epochs. Shown are power-law fits to the 
data: z/ a oc r°- 15±(U2 , Vja oc r 1 - 7 * ' 4 , and F m oc r a47±a11 . The fitted spectra are show in g, 
and the respective fit parameters can be found in Table 0. 



